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Background: It is now possible to map neural connections in vivo across the whole brain 
(i.e., the brain connectonne).This is a promising development in neuroscience since many 
health and disease processes are believed to arise from the architecture of neural networks. 

Objective: To describe the normal range of hemispheric asymmetry in structural connec- 
tivity in healthy older adults. 

Materials and IVIethods: We obtained high-resolution structural magnetic resonance 
images (MRI) from 17 healthy older adults. For each subject, the brain connectome was 
reconstructed by parcelating the probabilistic map of gray matter into anatomically defined 
regions of interested (ROIs). White matter fiber tractography was reconstructed from 
diffusion tensor imaging and streamlines connecting gray matter ROIs were computed. 
Asymmetry indices were calculated regarding ROI connectivity (representing the sum of 
connectivity weight of each cortical ROI) and for regional white matter links. All asymme- 
try measures were compared to a normal distribution with mean = 0 through one-sample 
Mests. 

Results: Leftward cortical ROI asymmetry was observed in medial temporal, dorsolateral 
frontal, and occipital regions. Rightward cortical ROI asymmetry was observed in middle 
temporal and orbito-frontal regions. Link-wise asymmetry revealed stronger connections 
in the left hemisphere between the medial temporal, anterior, and posterior peri-Sylvian 
and occipito-temporal regions. Rightward link asymmetry was observed in lateral temporal, 
parietal, and dorsolateral frontal connections. 

Conclusion: We postulate that asymmetry of specific connections may be related to 
functional hemispheric organization. This study may provide reference for future studies 
evaluating the architecture of the connectome in health and disease processes in older 
individuals. 
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INTRODUCTION 

With recent advances in structural neuroimaging, it is now pos- 
sible to track medium and large-scale pathways of white matter 
fibers and to construct the map of neural connectivity across 
the entire brain (the brain connectome) (1, 2). The brain con- 
nectome constitutes a promising development in neuropsychiatry 
since many physiological and pathological processes are believed 
to affect the architecture of neural networks (3-6). For example, 
normal cognitive development spanning from childhood to senior 
years is traditionally believed to be associated with maturation of 
neural networks supporting cognitive domains such as attention 
(7), memory (8), language (9), and executive function (10). Sim- 
ilarly, neurological diseases including epilepsy and dementia are 
associated with pathological rearrangements or impoverishment 
of normal networks at a systems level (11, 12). Furthermore, psy- 
chiatric diseases such as schizophrenia (13), bipolar disorder (14), 



and addiction (15) are related to reinforcements of pathological 
networks. 

In order to better understand how different biological processes 
affect the normal connectome, it is important to accurately charac- 
terize the connectome organization in healthy individuals. Specif- 
ically, it is important to define the degree of individual variability 
and hemispheric asymmetry that can be expected from the normal 
population. 

Since hemispheric asymmetry is abnormal in neuropsychiatric 
diseases such as schizophrenia (16), bipolar disorder (17, 18), and 
depression (19), in this study we aimed to describe the patterns 
of hemispheric asymmetry and individual variability of the brain 
connectome obtained from a cohort of healthy senior individuals. 
We employed high-resolution magnetic resonance images (MRI) 
to reconstruct structural brain connectivity based upon white 
matter pathways linking anatomically defined cortical regions of 
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interest (ROIs). We aimed to describe the group-wise distribution 
of asymmetries involving cortical connectivity and white matter 
pathways. 

MATERIALS AND METHODS 
SUBJECTS 

We studied 17 right-handed healthy subjects (mean age 53 years, 
SD = 7 years, range = 40-76) who were recruited from the local 
community. None of the subjects had a history of neurological, 
psychiatric, or chronic medical illnesses. All patients signed an 
informed consent to participate in this study. The Institutional 
Review Board at the University of South Carolina approved this 
study. 

IMAGE ACQUISITION 

AH subjects underwent MRI scanning at a 3T Siemens Trio 
equipped with a 12-channel head coil located at the Uni- 
versity of South Carolina, yielding: (1) Tl -weighed images 
(3D MP-RAGE, TR = 2250ms, TE = 4.15ms, 256 x 256 matrk, 
256 mm x 256 mm FOV, parallel imaging GRAPPA = 2, 80 ref- 
erence lines, TA = 377s); and (2) Diffusion weighted images 
(dMRI)-EPI scan (30-directions with fc= 1000s/mm2 and 
fo = 2000s/mm^, TR = 6100 ms, TE= 101ms, 82 x 82 matrix, 
222 mm X 222 mm FOV, parallel imaging GRAPPA = 2, 80, 45 
contiguous 2.7 mm axial slices, TA = 390 s). 

IMAGE PRE-PROCESSING 

The construction of the connectome involved two parallels pre- 
processing steps, namely, the segmentation of the cerebral cortex 
into multiple anatomical ROIs and reconstruction of white matter 
fibers. These steps are explained below: 

SEGMENTATION OF THE CEREBRAL CORTEX 

Tl -weighted MR images were converted into NIfTI format uti- 
lizing the dcm2nii tool from the MRIcron software package (20). 
Images in native space were non-linearly normalized into stan- 
dard MNI space using the Clinical Toolbox (21) employing unified 
segmentation-normalization routines as part of the software Sta- 
tistical Parametric Mapping (SPM8) . Of note, the Clinical Toolbox 
was particularly designed to accurately quantify tissue volumes in 
seniors and older adults (21). 

This step yielded probabilistic maps of gray and white matter 
in MNI space. 

Next, the base b = 0 T2 -weighted dMRI volume was linearly 
transformed to standard space utilizing a boundary-based regis- 
tration approach (22). This step was performed using FMRIB's 
Linear Image Registration Tool (FLIRT), as part of FMRIB Soft- 
ware Library (FSL) (23). The registration parameters were then 
used to transform the Automated Anatomical Labeling (AAL) 
atlas (24) and the white and gray matter probabilistic maps onto 
the dMRI space. Once in dMRI space, a map of cortical regions 
segmented according to AAL was obtained by overlaying the regis- 
tered AAL atlas onto the registered probabilistic gray matter map. 
The intersection between these images (including only voxels with 
a probability greater than 50% of being gray matter) represented 
the segmented cortical map. A list describing all ROIs used in this 
study can be observed in Table Al in Appendix. 



WHITE MATTER FIBER RECONSTRUCTION 

Extraction of diffusion gradients was performed with dcm2nii 
(20). The dMRI volumes were aligned to the b = 0 dMRI image 
using the FSL FLIRT tool (23). In diffusion space, whole brain trac- 
tography was reconstructed with the software Diffusion Toolkit 
(25) according to the following parameters: (1) angle thresh- 
old = 45°, (2) inclusion mask derived from the average of diffusion 
weighted signal and from the white matter probabilistic map reg- 
istered to dMRI space, (3) FACT propagation algorithm, (4) spline 
filter. 

CALCULATION OF THE CONNECTOME 

From each patient, the path of each tractography streamline was 
assessed. All streamlines were seeded in white matter, and the 
end-points of each streamline were computed. Streamlines with 
end-points within ROIs were counted as links between these ROIs. 
Streamlines with end-points outside ROIs were discarded. After all 
streamlines were assessed, the result was a weighted connectivity 
matrix A, where the entry Aij represented the number of stream- 
lines connecting regions and j (i.e., the weighted link between 
and j). Note that only direct links between regions and j were 
included in the link A^ . 

Finally, each link weight was corrected based on the surface 
of the connected ROIs and the distance between the ROIs, as 
proposed by Hagmann et al., where the link weight is inversely pro- 
portional to the sum of the linked ROI surfaces and fiber length, in 
order to account for tractography bias related to size of connected 
ROIs and distance traveled by the streamline (2). 

CONNECTOME ASYMMETRY 

Asymmetry was evaluated for ROIs and for links 

Regions of interested asymmetry was calculated by assessing the 
sum of link weights connecting an ROI, in comparison with the 
homologous ROI in the contralateral hemisphere. The connectiv- 
ity of each ROI was computed without discrimination regarding 
the opposite end of the streamline. For example, when computing 
the connectivity of the left hippocampus, all possible connections 
of the left hippocampus were computed, including connections 
to ROIs in the ipsUateral and contralateral hemispheres. Con- 
nections of the left hippocampus to itself or to the contralateral 
hippocampus were discarded. 

For each ROI, cortical connectivity asymmetry was calcu- 
lated according to the following asymmetry index (Al) = (R- 
L)/[(R-|- L)/2]; where L represents the connectivity of the ROI in 
the left hemisphere and R represents the connectivity of the ROI 
in the homologous ROI in the right hemisphere. A one-sample 
f-test was performed to evaluate whether the distribution of ROI 
asymmetries across all subjects was statistically different than a dis- 
tribution with mean = 0. This step was performed for each ROI. 
Note that ROI asymmetry represents the sum of connections to 
an ROI, irrespective to which other ROI is being linked to that 
ROI. Thus, ROI asymmetry should be interpreted in the context 
of link- wise asymmetry. 

Link-wise hemispheric asymmetry was calculated by assessing 
the difference in weight for each link, according to the same asym- 
metry index Al = (R-L)/[(R-|- L)/2]; where L represents the link 
between ROIs within the left hemisphere and R represents the 
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link between the homologous ROIs within the right hemisphere. 
Note that inter-hemispheric connections were excluded from this 
calculation. Cerebellar links were also excluded. Only supratento- 
rial links within the same hemisphere were assessed. A one-sample 
t-test was performed to evaluate whether the distribution of asym- 
metries for each link across subjects was statistically different than 
a distribution with mean = 0. 

RESULTS 

The average connectivity matrix is demonstrated in Figure 1. The 
ROIs are numbered in accordance with the glossary from Table Al 
in Appendix. Briefly, regions 1-45 represent ROIs located in the 
left hemisphere, whilst regions 46-90 represent ROIs located in the 
right hemisphere. As such, within the average connectivity matrix 
with 90 X 90 entries, the upper left quadrant demonstrates links 
between ROIs in the left hemisphere, and the lower right quadrant 
demonstrates links between ROIs within the right hemisphere. 
The right upper and the left lower quadrants of the connectivity 
matrix represent reciprocal left to right connections. Since con- 
nections are not directed (i.e., the strength of connectivity from 
region to j is the same as the connectivity from region jtoi), the 
matrix is symmetrical along its main diagonal. As expected, links 
within the same hemisphere exhibited a higher weight compared 
with connectivity links between different hemispheres. 



The average link- wise hemispheric asymmetry is also demon- 
strated in Figure 1. This matrix only includes 45 x 45 entries, 
where each cell represents the asymmetry index for each link. In 
order to avoid false estimations of asymmetry in links that were 
tracked only in a few subjects, we only included links that were 
tracked in greater than 75% of subjects (i.e., at least 13/17 sub- 
jects). Each matrix cell represents the asymmetry index for the 
connection between the ROI listed in the column and the ipsilat- 
eral ROI listed in the row. The ROIs are numbered in accordance 
with Table Al in Appendix (from 1 to 45, regardless of side). 

The distribution of ROI asymmetries is shown in Figure 2. The 
box-plot demonstrates the range of asymmetries across all sub- 
jects for each ROI. As expected, there were no regions of extreme 
asymmetry. A one-sample t-test revealed that among all possi- 
ble 45 ROIs, a significant leftward asymmetry was noted on the 
inferior occipital region, fusiform gyrus, lingual gyrus, parahip- 
pocampal gyrus, amygdala, inferior frontal operculum, superior 
frontal gyrus, and mid cingulate gyrus. Conversely, a rightward 
asymmetry was noted for the middle temporal gyrus and superior 
frontal orbital region. These results are shown in Table 1. 

Link-wise asymmetry was observed toward both hemispheres. 
A significant leftward asymmetry was noted on the following rec- 
iprocal connections: amydgala to parahippocampal gyrus; middle 
to inferior occipital regions; fusiform to lingual gyri; fusiform to 
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FIGURE 1 I Upper left panel: average connectivity matrix from all 
subjects - the scale bar represents ioglconnectivity weight). Upper 
right panel - link-wise average asymmetry index. Lower left panel - T 
values from a one-sample f-test evaluating the differences between the 
asymmetry indices and a distribution with mean = 0. Lower right 



panel - p values from the one-sample f-test. For all plots, the x- and 
y-axis represent ROIs numbered from 1 to 90 in accordance, where 
ROIs 1-45 are located in the left hemisphere, and 46-90 in the right 
hemisphere. The legend for ROI numbering can be seen InTableAl in 

Appendix. 
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FIGURE 2 I Box-plot demonstrating the distribution of asymmetry indices for all ROIs. Positive asymmetry indices represent a leftward asymmetry 



Table 1 | Asymmetry of cortical ROI global connectivity. 



Table 2 | Asymmetry of cortical links. 



Region Asymmetry index 


P 


T 


Link between 


Asymmetry index P 


T 




LEFTWARD ASYMMETR^^^^^^^^I 






ILeFTWARD ASYMMETRY \ 







Occipital inf 


-0.0381 


0.0023 


-3.6133 


Fusiforrm 


-0.0486 


0.0065 


-3.1304 


Lingual 


-0.0376 


0.0067 


-3.1154 


Parahippocampal 


-0.0320 


0.0073 


-3.0732 


Amygdala 


-0.0237 


0.0241 


-2.4916 


Frontal inf oper 


-0.0257 


0.0361 


-2.2873 


Frontal sup 


-0.0280 


0.0464 


-2.1587 


Cingulum mid 


-0.0232 


0.0488 


-2.1321 


RIGHTWARD ASYMMETRY ^^^^^^^^^^^^^^^ 


Temporal mid 


0.0368 


0.0100 


2.9201 


Frontal sup orb 


0.0445 


0.0283 


2.4112 



Only ROIs with a t-test p value less than 0.05 are shown. 

occipital inferior gyrus; insula to inferior frontal opercular region; 
fusiform to parahippocampal gyrus; precuneus to lingual gyrus; 
angular to superior parietal region; precuneus to mid cinculate 
gyrus. Conversely, a significant rightward asymmetry was noted 
on the following reciprocal connections: superior temporal gyrus 
to rolandic opercular region; middle temporal gyrus to middle 
occipital region; inferior to superior parietal regions; insula to 
inferior frontal triangularis region; precuneus to superior parietal 
region. These results are summarized in Table 2. 

These results are demonstrated in Figure 1, which shows the 
distribution of T scores and p values for all links, represented in 
45 X 45 matrices, where each matrix cell demonstrates the T value 



Amygdala 


Parahippocampal 


-0.8411 


0.0011 


-3.9764 


Occipital inf 


Occipital mid 


-0.8139 


0.0080 


-3.0298 


Fusiform 


Lingual 


-0.5496 


0.0087 


-2.9857 


Fusiform 


Occipital inf 


-0.7223 


0.0113 


-2.8617 


Insula 


Frontal inf oper 


-0.7726 


0.0156 


-2.7069 


Fusiform 


Parahippocampal 


-0.4685 


0.0239 


-2.4952 


Precuneus 


Lingual 


-0.6617 


0.0406 


-2.2281 


Angular 


Parietal sup 


-0.6378 


0.0421 


-2.2092 


Precuneus 


Cingulum mid 


-0.5340 


0.0453 


-2.1715 


RIGHTWARD ASYMMETRY ^^^M 


Temporal sup 


Rolandic oper 


0.6482 


0.0040 


3.3592 


Temporal mid 


Occipital mid 


0.6424 


0.0098 


2.9310 


Parietal inf 


Parietal sup 


0.6973 


0.0113 


2.8605 


Insula 


Frontal inf tri 


0.4912 


0.0129 


2.7959 


Precuneus 


Parietal sup 


0.6894 


0.0157 


2.7016 



Only links with a t-test p value less than 0.05 are shown. 

or the p value for the asymmetry distribution for the connection 
between the row ROI and the column ROI. Figure 3 demon- 
strates the anatomical distribution of links with an absolute mean 
asymmetry higher than 0.5. The three-dimensional anatomical 
reconstruction of regional links was defined based on an in-house 
developed atlas of anatomical connectivity involving all 90 ROIs 
used in this study. The location of travel of all streamlines connect- 
ing each possible pair of ROIs was defined based on the spatial 
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FIGURE 3 I Anatomical demonstration of the location of linl<s witli an absolute asymmetry index greater than 0.5 Links in red exhibit a leftward 
asymmetry, while links in blue demonstrate a rightward asymmetry. 



distribution of center points (or "centroid") of serial transverse 
sections across the white matter streamlines corresponding to 
each link, through an in-house modified version of the meth- 
ods described by Garyfallidis et al (26). The centroids for each 
link were connected to define the main pathway of streamline 
travel. This step was repeated for all possible links. The utility of 
this approach was exclusively for anatomical visualization of the 
results, as demonstrated in Figure 3. 

DISCUSSION 

In this study, we investigated the individual variability of hemi- 
spheric asymmetry in structural connectivity in healthy senior 
individuals. By reconstructing the structural connectome from 
each subject, and by assessing the distribution of asymmetry 
indices related to cortical ROI connectivity and link weight, 
we observed that, within a sample composed of healthy older 
adults, there was mild but noticeable hemispheric asymmetry in 
structural connectivity. 

We observed a more prominent leftward asymmetry in cor- 
tical connectivity, i.e., a larger number of regions demonstrated 
a higher degree of connectivity in the left hemisphere. Specifi- 
cally, cortical ROIs located in the occipital lobe, medial tempo- 
rsJ, dorsolateral and medial frontal lobe and cingulate exhibited 
a higher weight of connectivity in the left hemisphere. Con- 
versely, fewer cortical regions demonstrated a rightward asymme- 
try, with the middle temporal gyrus and the orbito-frontal regions 
exhibiting a significantly higher cortical connectivity on the right 
hemisphere. 

In turn, we also observed regional hemispheric asymmetry in 
relationship with the strength of connectivity between specific 
ROIs. In accordance with the previous observation about global 
cortical connectivity, a leftward asymmetry was also more com- 
monly observed among links. Interestingly, a leftward asymmetry 
was observed on peri- Sylvian and medial temporal - occipi- 
tal regions. Conversely, a rightward asymmetry was noted on 
lateral temporal - frontal - occipital and parietal links. While 
did not test the relationship between these links and cogni- 
tive performance, it is possible to speculate that some lateral 
asymmetry may be related to functional specialization of some 
of these connections. For instance, leftward asymmetry may 
be associated with dominant hemisphere language processing 
involving verbal memory (parahippocampal - amydala connec- 
tions) (27), phonological processing (insula - frontal opercu- 
lum connections) (28), and semantic retrieval (29) (angular 



gyrus - superior parietal connections). Conversely, a rightward 
asymmetry may be observed in networks associated with visual- 
spatial processing (30) (intraparietal connections, precuneus- 
parietal connections). 

It should be noted that we adopted a liberal statistical thresh- 
old (i.e., the level of statistical significance from the one-sample 
f-tests was not corrected based on multiple comparisons). We 
adopted a liberal threshold since we expected that the degree of 
asymmetry exhibited by our population would be mild, and given 
the number of multiple comparisons, a more stringent thresh- 
old would preclude the evaluation of the locations with a higher 
degree of asymmetry. Nonetheless, given our sample size, it is 
possible that some of our observed asymmetries may be related 
to sample bias and may constitute false positives. For this rea- 
son, we recommend the interpretation of the results from this 
manuscript in this context. The asymmetry index is possibly 
a better representation of the magnitude of asymmetry, rather 
than over- emphasizing the importance of links or ROIs with 
p < 0.05. 

The results reported in this study should help the contextual 
evaluation of other connectome studies applied to health and dis- 
ease. We demonstrated the range of asymmetry in a small cohort 
of normal older adults with the purpose of providing a reference 
for future studies evaluating processes that affect neural network 
organization. Thus, future studies should also be interpreted with 
special attention to specific characteristics of the demographics 
from the population studied. 
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APPENDIX 



Table A1 | Legends for ROI numbering. 



Number 


Side 


ROI 


Number 


Side 


ROI 


1 


Left 


Precentral 


46 


Right 


Precentral 


2 


Left 


Frorntal sup 


47 


Right 


Frontal sup 


3 


Left 


Frontal sup orb 


48 


Right 


Frontal sup orb 


4 


Left 


Frontal mid 


49 


Right 


Frontal mid 


5 


Left 


Frontal mid orb 


50 


Right 


Frontal mid orb 


6 


Left 


Frontal inf oper 


51 


Right 


Frontal inf oper 


7 


Left 


Frontal inf tri 


52 


Right 


Frontal inf tri 


8 


Left 


Frontal inf orb 


53 


Right 


Frontal inf orb 


9 


Left 


Rolandic oper 


54 


Right 


Rolandic oper 


10 


Left 


Supp motor area 


55 


Right 


Supp motor area 


11 


Left 


Olfactory 


56 


Right 


Olfactory 


12 


Left 


Frontal sup medial 


57 


Right 


Frontal sup medial 


13 


Left 


Frontal mid orb 


58 


Right 


Frontal mid orb 


14 


Left 


Rectus 


59 


Right 


Rectus 


15 


Left 


Insula 


60 


Right 


Insula 


16 


Left 


Cingulum ant 


61 


Right 


Cingulum ant 


17 


Left 


Cingulum mid 


62 


Right 


Cingulum mid 


18 


Left 


Cingulum post 


63 


Right 


Cingulum post 


19 


Left 


Hippocampus 


64 


Right 


Hippocampus 


20 


Left 


Parahippocampal 


65 


Right 


Parahippocampal 


21 


Left 


Amygdala 


66 


Right 


Amygdala 


22 


Left 


Calcarine 


67 


Right 


Calcarine 


23 


Left 


Cuneus 


68 


Right 


Cuneus 


24 


Left 


Lingual 


69 


Right 


Lingual 


25 


Left 


Occipital sup 


70 


Right 


Occipital sup 


26 


Left 


Occipital mid 


71 


Right 


Occipital mid 


27 


Left 


Occipital inf 


72 


Right 


Occipital Inf 


28 


Left 


Fusiform 


73 


Right 


Fusiform 


29 


Left 


Postcentral 


74 


Right 


Postcentral 


30 


Left 


Parietal sup 


75 


Right 


Parietal sup 


31 


Left 


Parietal inf 


76 


Right 


Parietal inf 


32 


Left 


Supramarginal 


77 


Right 


Supramarginal 


33 


Left 


Angular 


78 


Right 


Angular 


34 


Left 


Precuneus 


79 


Right 


Precuneus 


35 


Left 


Paracentral lobule 


80 


Right 


Paracentral lobule 


36 


Left 


Caudate 


81 


Right 


Caudate 


37 


Left 


Putamen 


82 


Right 


Putamen 


38 


Left 


Pallidum 


83 


Right 


Pallidum 


39 


Left 


Thalamus 


84 


Right 


Thalamus 


40 


Left 


HeschI 


85 


Right 


HeschI 


41 


Left 


Temporal sup 


86 


Right 


Temporal sup 


42 


Left 


Temporal pole sup 


87 


Right 


Temporal pole sup 


43 


Left 


Temporal mid 


88 


Right 


Temporal mid 


44 


Left 


Temporal pole mid 


89 


Right 


Temporal pole mid 


45 


Left 


Temporal inf 


90 


Right 


Temporal inf 
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